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Kinetics of merging and breaking down vortex loops is the important part of the whole vortex 
tangle dynamics. Another part is the motion of individual lines, which obeys the Biot-Savart law 
in presence of friction force and of applied external velocity fields if any. In the present work we 

■ evaluate the coefficients of the reconnection rates A(h, h, I) and B(l, h, h). Quantity A is a number 
(per unit of time and per unit of volume) of events, when two loops with lengths iiand I2 collide 
and form the single loop of length I = h + h. Quantity B(l,li,h) describes the rate of events, 

00 ■ when the single loop of the length I breaks down into two the daughter loops of lengths h and h- 

These quantities ave evaluated as the averaged numbers of zeroes of vector S s (£2, £1, i) connecting 
two points on the loops of £2 and fi at moment of time t. Statistics of the individual loops is taken 
from the Gaussian model of vortex tangle. 
PACS-number 67.40 

O 
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■ I. INTRODUCTION 
Of 

Quantized vortices appeared in quantum fluids and other systems play a fundamental role in the properties of the 
latter. For that reason they have been an object of intensive study for many years (for review and bibliography see e. 
g. [1]). The greatest success in investigations of dynamics of quantized vortices has been achieved in relatively simple 
cases such as a vortex array in rotating helium or vortex rings. However these simple cases are rather exception than 
1 . ■ a rule. Due to extremely involved dynamics initially straight lines or rings evolve to form highly chaotic structure, 
so called vortex tangle (see e.g. papers [2], [3] and review [4]). The vortex tangle consists of a set of vortex loops 
q . of different lengths. The individual vortex loops evolves obeying the Biot-Savart law, undergoing the friction force 
and the applied external velocity field. All these processes are good understood, although the according equations of 
*~ ~ 1 motion arc highly nonlinear (and nonlocal) therefore general analysis can be made with use of numerical methods. 

Besides motion of each of the individual loops there is one more element of general evolution of the vortex tangle 
related to collision of loops, or intersection of elements of vortex lines. During intersection of lines the very complicated 
f— ( ■ process, related to arrangement of the vortex core takes place [5] . Nevertheless this process is relatively short, therefore 
it is usually accepted that the filaments instantly reconnect whenever they intersect each other. It is widely appreciated 
that the reconnection processes influence both the structure and evolution of vortex tangle. However the questions 
how it happens and what mechanisms are responsible for these remain open. For instance, Feynman in his pioneering 
paper devoted to superfiuid turbulence proposed scenario how the vortex tangle decays in absence of applied external 
counterflow. According to this scenario a fusion of small vortex rings into larger ones as well as a breakdown into 
smaller ones is possible at the moments of the reconnection events (see Fig. 1). In assumption that on the average 
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FIG. 1: The reconnection process schematically (Feynman, 1955, Fig. 10). (a) initial stage, (b) and (c) are stages of collapse, 
(d) reconnection stage, (e) stage of degeneration of vortex rings into thermal excitations. 
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FIG. 2: Schematic sketch of the fusion and breaking-down of vortex loops. Rates of these processes characterized by the 
coefficients A(h,h,l) and B(l,h,h). 

the last property dominates i.e. the cascade like process of formation of smaller and smaller loops forms. When the 
scale of the small rings becomes of the order of the interatomic distances, which is the final stage of the cascade, the 
vortex motion is degenerated into thermal excitations. This key idea that degeneration of the vortex tangle occurs 
due to cascade-like transferring of the length in space of scale of sizes of vortex loops was indirectly confirmed only 
in numerical calculations, where the procedure of artificial elimination of small loops had been used [6]- [9]. 

In spite of the recognized importance of the reconnecting loops kinetics, the numerical results remain main source 
of information about this process. The obvious lack of theoretical investigations interferes with deep insight in the 
nature of this phenomena (this question had been recently discussed in [10]). For instance it is not clear how the 
cascade of length in space of vortex loops sizes is formed, what mechanisms are responsible for this, what quantities 
determine an intensity of cascade, and why at all breakdown of the loop prevails. 

Of course the scarcity of analytic investigations related to incredible complexity of the problem. Indeed we have 
to deal the set of objects with not fixed number of elements, they can born and die. Thus, some analog of the 
secondary quantization method is required with the difference that objects (vortex loops) themselves possess an 
infinite number of degree of freedom with very involved dynamics. Clearly this problem can hardly be resolved in 
nearest future. Recently in [11] much more modest approach, based on the "rate equation" for distribution function 
n(l) was elaborated in context of cosmic strings. Following this work we introduce distribution function n(l,t) of 
density of loop in "space" of their lengths. It is defined as a number of loops (per unit of volume) with lengths 
between I and I + dl. Due to reconnection processes n(l,t) can vary. 

We discriminate two types of processes, namely the fusion of two loops into the larger single loop and the breakdown 
of single loop into two daughter loops. The kinetic of vortex tangle is affected by the intensity of the introduced 
processes. The intensity of the first process is characterized by the rate of collision A(h,l2,l) of two loops with 
lengths Ziand h and forming the loop of length I = l\ + 1 2 . The intensity of the second process is characterized by the 
rate of self-intersection B(l, I1J2) of loop of the length / into two daughter loops with lengths l\ and I2 (see Fig. 2). In 
view of exposed above we can directly write out the master " kinetic" equation for rate of change distribution function 
n(l,t). 



All of the processes are depicted at the left of each line. Clear that the "kinetic" equation has "bookkeeping" character. 
Physics of this approach lies in the "correct" determinations of coefficient of that equations A(li,l 2 J) and B(l, i 1} i 2 ) 
on the base of some more or less plausible model. 



dn(l,t) 
dt 



(1) 




FIG. 3: Schematic sketch of vortex loop. Elements of line are described as vectors s(£), where the label variable £ is taken here 
as the arc length. We associate the moment of intersection with the vanishing of vector S(£i,£i,t) connecting points s(^2,t) 
and s(£i, t). 

In the present work we derive some general (independent on the model) relations for reconnection rates A(l\, I2, 1) and 
B(l, li, I2)) ( Section II). Then in Sec. Ill we apply these relations to evaluate the reconnection rates on the base of 
Gaussian model of the vortex tangle elaborated earlier [12]. Section IV is devoted to conclusion and plans for future 
investigations. 



In this section we will formulate mathematical definition for quantities A(Ji,l2,l) and B(l, l\,l2)- We will start 
with the latter quantity B(l,li,l 2 ). By definition its physical meaning is just frequency of events when part of line 
with total length I intersects the same line to create two two daughter loops with lengths l\ and I2 (see Fig. 2), that 
is self-crossing event. As it is already stated we assume that each crossing event leads to reconnection of lines. The 
elements of vortex line are described as a function s(£, t) (See Fig. 3) that is time dependant radius-vector of the 
points resting on loop. Variable £ labels the points of the loop. It is convenient to choose variable £ to be equal to 
the arc length, (0 < £ < I). Let us consider function 



that is just 3D vector connecting points s(^ 2 ,t 2 ) and s(£i, ti). Clearly that condition S;,(£2, £1, t) = implies that the 
self-crossing event of parts of line with label-coordinates £2, £1 occurs at moment of time t. The quantity S{,(£ 2 , £1, i) 
is fluctuating 3-component function of three arguments £2 ; £1 ? ^- We are interested in how often Sb(£2,£i,i) vanishes 
in cube of space £ = {£2, £1, t}. Points in space £ = {£2, £i 5 1} where function St>(£2, £1, t) vanishes, are points of cross 
of 3 surfaces Si (£2, £1, t) = 0, i = x, y, z in space £ = {£2, £1, t} as it is shown in Fig. 4. 

From theory of generalized function it follows that number of these points (we denote them below as £ a ) can be 
expressed via S- function of quantity Sf,(£2,£i,t) with the help of following formula . 



Here X, Y, Z are the components of vector S{,(£2, £2, £1, 0) . By integration of both parts of (3) over d£ic?£2 we would 
obtain the full number of intersections (per unit time). If, further to introduces additional constraint ^(£2 — £1 — h) 
implying that the distance (along line) between chosen points is equal to I and to integrate over d£id£2, we obtain 
the rate of self-intersection of line with length I and breakdown it into pieces l\ and I — l\. 



II. MATHEMATICAL IDENTITIES FOR A(h,h, I) AND B(l, h, l 2 ). 



S s (£ 2 ,£i,i) =s(£ 2 ,i) 



s(£i,i), 



(2) 




(3) 



FIG. 4: The reconnection events can be described as zeroes of function St(^2, £1, t). In space of its variables £ = {£2, £1, i} these 
points are cross of 3 surfaces <Si(£2, £1, i) =0, i = x,y,z. 



In addition we have to do averaging over all possible fluctuating configurations. Thus the coefficient B(l\, I — h, I) 
with dimension [s] — s _1 cm _1 is equal to 

B(h,l-h,l) = J J deid&*(6 - £ x - Zi) 

To obtain coefficient A(l\,l2,l) we use the similar procedure. Let us consider two loops with length l\ and I2. Our 
purpose now to find the rate A(l\, I2, 1) of fusion of these two loops into one loop of length I = li + 1%. Dimension of 
A{li,l 2 ,l) is [A] — cm 3 s~ 1 . As previously we describe vortex filaments by positions of radius vectors their elements 
s(£i, t) and s(£2, t). Here we have two label variables £1, £2 belonging to different loops and running in limits (0 < £ < 
li) and (0 < £ < 1 2) respectively. One more important difference with the previous case is that both functions s(£ 1; i) 
and s(^2,t) should depend on "initial" positions s(£2 = 0, t) — Ri(t) and s(^2 = 0, t) = R2(t), chosen arbitrary. Of 
course in previous case of self-intersection of single loop, quantity s(£,i) also depended on "initial" positions R(t), 
but it did not influence the rate of self-intersection. Now for case of the fusion this dependance is important, since 
very distant loops have small probability to collide. Let us introduce the "fusion" functions 

S / (6,a,ii)=s(6,t)-s(£i,t). (5) 

Repeating the considerations for case of the single loop we find that number of reconnection (per unit of time) of 
points £2,^1 formally coincides with (3) 



d(X,Y,Z) 
5(£ 2 ,£i,t) 



*(S 6 (6,£i,t)) 



(4) 



c=c„ 



with the difference that £2,£i belong to different curves. Since intersections of any elements of lines lead to fusion of 
the loops we have to integrate (6) over tf£id£2. The result obtained is valid for chosen pair of loops. To obtain the 
total number of events we have to multiply the result obtained by quantity n(li)n(l 2 )dllidll2, which is just a full 
number of loops of chosen sizes in the whole volume. Comparing with the master kinetic equation (1) we find the 
final expression for fusion coefficient A(l\, I2, 1) 



d(X,Y,Z) 
d(£ 2 ,£i,t) 



*(S/(&,&,*)) 



(6) 



c=e„ 



Thus we obtained expressions (4) and (7), which allow to calculate the rates of reconnections of the fusion and 
breakdown vortex loops. They are however just formal mathematical identities. Concrete results depend on statistics 
of individual lines. Therefore to move further we have to ascertain the procedure for averaging. It will be done in the 
next section. 
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FIG. 5: Snapshot of the "average" vortex loop obtained from analysis of the statistical properties. Close (A£ <JC R) parts of 
the line are separated in 3D space by distance A£. The distant parts (7? <C A£) are separated in 3D space by the distance 
\/£<>A£, i.e. the vortex loop has the typical random walking structure. 



III. GAUSSIAN MODEL CASE 



A. The generalized Wiener distribution 



To evaluate quantities B{1\, I2, 1) and A{1\, I2, 1) written in form (4) and (7) one needs to know statistics of individual 
loops. In general case this statistics should be extracted from investigation of the full dynamical problem. The 
according statement of such problem includes equation of motion (Biot-Savart law for quantum vortices), dissipative 
effects (interaction with normal component) and additional Langevin force responsible for chaotic behavior. The 
problem becomes more involved by circumstance that other loops also influence dynamics of the chosen loop. At 
this stage we choose another way, namely we use the Gaussian model of the vortex tangle elaborated by author [12]. 
To develop this approach the trial distribution functional of Gaussian form had been consructed. This functional 
absorbed all properties of the supcrfluid turbulence known from both experimental studies and numerical simulations. 
According to this model the "average" vortex loop has a typical structure shown in Fig. 5. The close parts of the 
loop separated (along line) by distance £2 — £1 smaller then the mean radius of curvature £0 are strongly correlated, 
(s'(£i, i)s'(£ 2 , t)) — > 1, (s' is the tangent vector) and line is smooth. Remote parts of the line (£2 — £,1 ~> £0) are 
not correlated at all, (s'(£i, i)s'(£ 2 , t)) — > 0. Thus for large separations the vortex loop has a typical "random walking 
structure" with the Wiener distribution. This " semifractal" behavior satisfies to the generalized Wiener distribution. 
Namely, the probability V({s(£,t)}) to find some particular configuration {s(£, t)} is expressed by the probability 
distribution functional (see for details in [12]) 



Here Af is normalizing factor, I is the length of curve, it is supposed to be much larger then mean radius of curvature 
£o- The probability distribution functional (8) introduced in [12] distribution considered anisotropy and polarization 
of vortex tangle in counterflowing turbulent helium. This fact is reflected by the circumstance that A"^(£i — £2) is 
the full matrix. In practice it is more convenient to deal with characteristic functional M / ({P(C, t)}) defined as 



P({s(£,t)})=A/-exp 




s /a (£i,i)A Q/3 (6 " 6)3^2,^1^2 



(8) 





The characteristic functional enables us calculate any averages depending on vortex lines configuration {s(£, t)} by 
simple functional differentiation. For instance the average tangent vector (s^(£i)) or the correlation function between 
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orientation of the different elements of the vortex filaments ^s^,(£i)s^(£ 2 )y are readily expressed via characteristic 
functional accordingly to the following rules: 



Km 



sw 



i6P a (tii) 



p = o 



^(6^(6)) 



S 2 W 



(10) 



p = o 



Calculation of the characteristic functional W({P(£, i)}) (9) on the base of the probability functional (8) is reduced 
to functional integration, which, in turn, reduces to the " full square procedure" . The result is 



l i 



(ii) 







To avoid unnecessary lengthy calculations we simplify the model expressed by the probability distribution functional 
(8), namely we omit both the anisotropy and polarization. In this case the matrix N af) (^ - £ 2 ) used in ([12]) can 
be taken as 



1 



6(1-2WW0 



exp 



(6-6 



26 

l 



(12) 



For small separation (£2 — 6 "C 6) i J2 a ^"(6 _ 6) — * 1/2 that guarantees that (s'(£, i)s'(£, t)) = 1, as it should be 
for smooth lines (we recall that for pure random walking line the quantity (s'(£, t)s'(£, t)) does not exist at all). For 
large separation (£ 2 — 6 ^ 6) exponents tends to (5 ce( 3(2y / 7r^o/6)(5(^i — £2), and correlation between tangent vectors 
weakens, (s'(£, t)s'(£, t)) — > 0. The second term in parenthesis of expression (12) appears due to closeness of the loops. 
It ensures that 

l 1 



J J d6<*6 (s'(6,t)s'(6,t)) =((sG,t)-s(0,i)) 2 ) =0 





as it should be for the closed lines (see explanations in ([12])). Thus the characteristic functional with function N 
satisfies to necessary " scmifractal" behavior of line and will be used further for evaluation of the rate reconnection 
coefficients. Here it is necessary to do two important comments. First one concerns the " semifractal" structure of 
the vortex loop. In principle, reconnections events occur due random walking behavior with non-zero probability 
of the meeting of any two elements. However, to evaluate the according rates on base of relations (4), (7), we 
need to know not only an average of <5(S(£2, 6> ^)) but also an average of the Jacobian. But the latter includes 
averaged derivatives (with respect two both label variable^, 6 and time). It is known fact however that quantities 
like (|s^(6)|) , ^ s L(6) s /j(6)^ e tc- are ill-defined for pure random walking structure, For that reason we need to 
smooth away the curve on small scales. The second remark concerns the mean radius of curvature £o- This quantity 
appears here as important parameter of theory of random walk. It is frequently called as " an elementary step" or 
persistence length. In fact the theory described by the Wiener distribution is failed for small scales less then £0 , 
therefore usually this value appears as a low cut-off of the whole approach. 



IV. EVALUATION OF A(h,l 2 ,l) AND B(l,h,h). 



A. evaluation of B(l,li,h) 



Let us come back to the rate coefficients. We start with the self-intersection processes. Positions of line elements 
s(£2)£),s(6)£) an d relative vector Sb(^2,£,i,t) are strongly fluctuating quantities having the Gaussian statistics. Due 
to the Wick theorem the average in integrand of (4) can be taken as a sum of all possible pairs of quantity 
S&(6>*>6>*) and its derivatives. Because of uniformity in £ space, quantity Sfc(^2,6^) depends on [£2 — 61 i f° r 
that reason all averages of structure ((dX/d^i)5(X(^2, 6> *))) vanish, therefore only pairs separately from S(,(^2, 6) *) 
and from its derivatives survive. As a result the average of production is equal to production of averages and each 
of the factors can be evaluated separately 

W(S 6 (6,6,i))>- (13) 



d(X,Y,Z) 

9(6, 6, t) 



<5(s fe (6,6,*)) 



C = Ca 



d(X,Y,Z) 

9(6,6,*) 
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As mentioned, the use of the characteristics functional (9), (11) allows to calculate any averaged functional of config- 
urations {s(£,i)} . Let us show how to evaluate (#S S (£2, £i, t)). With use of the standard integral representation for 
(5-function 

s ^=(k £«-*• 

we rewrite (5Sb{&,t, £i,i)) as 



(<5S 6 (£ 2 ,£i,i)) = ^3 J (cxp[*y(s(C 2 ,t 2 )-s(6,ii))])d 3 y 



(2tt) 3 



exp \i / ys'(£,i)d£<ft )d 3 y. 



(14) 



Comparing (9) and (14) we conclude that the integrand in last term of (14) is just the characteristic functional 
W({P(£,t)}), taken at value of P(£,t) 



p(£) = -y*(£-£i)0(6-O- 



(15) 



Here 0(£) is the unit step-wise function. Relation (15) implies that we choose in integrand of the characteristic 
functional only points lying in interval from £i to £2 on the curve. Substituting this value P(£, i) expressed by (15) 
in CF (11) with function iV Q/3 (£i - £2) in form (12) we obtain 



(£S 6 (£ 2 ,i,£i,i)) = 



(2tt) 2 



exp 



-y 2 ^ ((6 



a)- -6)70 



d 3 y = 



27r3/%((£ 2 -a)-(6-6) 2 /0 



3/2 



( 16 ) 

Evaluation of absolute value of Jacobian in (13) we perform by use of | J| = \fiP. Furthermore it would be convenient 
for the sake of generalization to use vector velocity V; = ds/dt, instead of calculation of it in explicit form, expressing 
velocity via the vortex filament configuration {s(£)}. Calculation of J 2 can be fulfilled writing Jacobian in explicit 
form and subsequent applying the Wick theorem. Simple but tedious calculations leads to result that 



J 2 = 2 (Vl) ((ds y /d^f) ((ds z /d&) 2 ) + 2 (Vl) ((ds v /d&) 2 ) ((dsz/dbf) + p.p. 



(17) 



where p.p. all permutations with respect to x, y, z. Estimating ((<9s y /<9£i) 2 ) and similar terms as 1/3 we obtain that 
I J I = I |V/| . After use of the integration J J d^id^S^ — £1 — ^1) (see 4) we finally obtain 



Bih,l-lul) = b.*V< m J m)i/a (18) 

where constant b s — (V3/8tt 9 / 4 ) w 0.0164772. We introduced in coefficient B the additional factor 1/2 in order 
to avoid the over-counting of the reconnection events, since decays I — > l\ + I2 and I — ► i 2 + describe the same 
process, though the both enter into equations. Let us recall that the quantity £ is the mean radius of curvature. 
The numerical factor b s ~ 0.0164772 is the result of particular approximation used in paper ([12]). If for instance we 
used pure Wiener distribution with the persistency length £0 and disregard the closeness of line, we would have for 
coefficient B{l\,l — h,l) the expression 

Bc(h,l-h,l)=b s *V l 7 - 1 (19) 

where b s 0.11. This result is very similar to result obtained earlier in paper [11] from the qualitative consideration. 
Yet it is important also that in this way it is not clear how to relate the persistency length £ with the mean radius 
of curvature. 
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B. evaluation of A(h, h, I) 



Let us now evaluate quantity A(l\, l 2 , 1) denned by relation (7). We again (as for the previous case) evaluate average 
from Jacobian and 5- function separately. Contribution from Jacobian coincides with the previous result \J\ = | |Vj|. 
The rest (5-function part can be evaluated with the help of the CF obtained above. Unlike previous case we have to 
know two-loop distribution function. Since we omit interaction of loops (until the reconnection event occurs) the CF 
for two loops with lengths Ziand l 2 is just production of the expressions of type (11) 

-J J p^aX'ta - 6)p^(6)^i^ 2 1 x (20) 







Quantities A r " /3 (£i — £ 2 ) an d N^fa — £ 2 ) differ from each other only by lengths of loops ^and l 2 , entering expressions 
for N" p . Further, by use of the standard integral representation for 5-function we have 

(5S f (t2,ti,t)) = ^3 J (cxp^y(s 2 (6,t 2 )-s 1 (ei,ii))])rf 3 y. (21) 

We stress again that the label variables £ 2 and £1 belongs two different loops. Let us introduce initial points si(0) 
and s 2 (0) and rewrite (21) in the following form: 

^3 J exp [-iy(s 2 (0) - Sl (0))] (exp [iy(s 2 (&) - s 2 (0))] exp [-iy( Sl (6) - Sl (0))]) d 3 y. 

Identifying further the "initial" positions si(0),s 2 (0) with quantities Ri,R 2 in formula (7) we rewrite it as 

A{lM) = VOdilME^ 1 / / rfRlrfR2 / / ^ 2 (22) 
J exp [-iy(Ra - Ri)] (exp [iy(s 2 (&) - s 2 (0))] exp [-iyfafa) - s^O))]) d 3 y. 

Let us introduce variables Ri — R 2 , (Ri + R 2 )/2. Integration over R 2 Ri gives 5(y), integration over (Ri + R 2 )/2 
gives the total volume of system. Further, integration over y gives unity, and integration over £1, £ 2 gives the production 
?iZ 2 . Thus we obtain the remarkable result, that for noninteracting loops quantity A(l±, i 2 , 1) responsible for merging 
of loops does not depend on statistics of individual loop at all and is equal 

A(l 1 ,l 2 ,l) = b m V l l 1 l 2 (23) 

Here b rn — 1/3. As earlier we introduced additional factor 1/2 to avoid the over-counting of the reconnection events. 

Results (23) and (18) (with not well determined factors b s and b m ) were also obtained in paper [11]. Authors used 
some qualitative picture of moving and colliding elements of lines. This fact confirms the validity of approach made 
in our work, which allows to use it for more complicated (in comparison with random walk) cases. 



V. CONCLUSION 



In the paper we discussed a role of kinetics of the merging and breaking down vortex loops in the whole vortex 
tangle dynamic. Following the work [11] we introduced the reconnection rates A(li, l 2 , 1) and B(l, l\,l 2 ). Quantity A 
is a number (per unit of time and per unit of volume) of the fusion of two loops of lengths li&nd l 2 with formation 
of the single loop of length I = l\ + l 2 . Quantity B(l,l\,l 2 ) describes rate of the breakdown of the loop with the 
length I into two daughter loops with lengths l\ and l 2 . We developed mathematical formalism, which enables us to 
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calculate A(li,l2,l) and B(l,li,l2). Briefly, these quantities can be found as an average number of zeroes of vector 
SSfefej £ l5 1) connecting points of £2 and £1 at moment of time t. The averaging procedure should be fulfilled by using 
statistics of individual loops, expressed e.g. by the probability distribution functional V ({s(£ , t)}) . Wc applied this 
formalism to evaluate the reconnection rates A(li,l2, and B(l, h,h) in case when vortex loop is semifractal object, 
which is smooth for small separations along the line and has a random work structure for distant (along line) points. 
Relations (7) and (4)are the key results of the paper. They are formal, and do not depend on statistic of loop. They 
can be applied for any physical situation provided that we are able to get statistics from the dynamical problem. In 
that sense the application of the results (7) and (4) for Gaussian case is just illustration how to use it in concrete 
sense. One of serious problem in this illustration is that in Gaussian approach the low scale was restricted by mean 
radius of curvature. In real vortex tangle £ m superfluid turbulent Hell this scale exceeds greatly the scales at 
which the real dissipativc mechanisms occur. Thus, the study of what would happen for scales smaller of mean 
radius of curvature is on of actual problems. Nevertheless results expressed by relations (23) and (18) are significant, 
they allows to analyze the kinetic equation (1), the attempts were made in papers [11], [13]. We also will apply these 
result in the forthcoming paper for study of the kinetic equation. 

I am grateful to participants of the workshop " Superfluidity under Rotation" (Manchester, 2005) for useful discussion 
of the results exposed above. This work was partially supported by grant N 03-02-16179 of the Russian Foundation 
of Fundamental Research. 
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